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Abstract 



The Discrete Nonlinear Schrodinger Equation is used to study the formation 

> ' 

lO ' of stationary localized states due to a single nonlinear impurity in a Caley 

(N 

' tree and a dimeric nonlinear impurity in the one dimensional system. The 

o , , 

\^ ' rotational nonlinear impurity and the impurity of the form — X I C T where a 

is arbitrary and x is the nonlinearity parameter are considered. Furthermore, 
S ' \ C \ represents the absolute value of the amplitude. Altogether four cases 

pH ■ are studies. The usual Greens function approach and the ansatz approach 

o ■ 

are coherently blended to obtain phase diagrams showing regions of different 

V~j ■ number of states in the parameter space. Equations of critical lines separating 

■ various regions in phase diagrams are derived analytically. For the dimeric 

problem with the impurity — x | C 1°^, three values of | Xcr \, namely, | 

Xcr 1= 2, at fj = and | Xcr |= 1 aud | for a = 2 are obtained. Last 

two values are lower than the existing values. Energy of the states as a 

function of parameters is also obtained. A model derivation for the impurities 

is presented. The implication of our results in relation to disordered systems 

comprising of nonlinear impurities and perfect sites is discussed. 
PACS numbers : 71.55.-i, 72.10.Fk 
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I. INTRODUCTION 



The Discrete Nonlinear Schrodinger equation (DNLS) is a set of N coupled differential 
equations. 

. dCm 



''X.mfm{\ Cm \)Cm ~\~ '^m,m-\-lC m+1 ~\~ ^m,m—lCi 



m— 1 



dt 

where Vm,m+i = Ki^+i,™; and m = 1, 2, 3, A^. (1) 

In eq.(l) the nonlinearity appears through functions fm{\ Cm |) and Xm is the nonlinearity 
parameter associated with the m-th grid point. Since, J2m I Cm P is made to unity by 
choosing appropriate initial conditions, | Cm P can be interpreted as a probability of finding 
particle at the m-th grid point. One way to derive this set of equations is to couple in the 
adiabatic approximation the vibration of masses at the lattice points of a lattice of N sites 
to the motion of a quasi particle in the same lattice. The motion of the quasi particle is 
described, however, in the frame work of a tight binding Hamiltonian (TBH). Same type 
of equation can also be obtained by nonlinear coupling of anharmonic oscillators through 
both positions and momenta of the oscillators p|-p!3|. The set of equations, thus derived, 
are called the discrete self-trapping equation (DST). These equations also posses a constant 
of motion analogous to X^m I Cm P in the DNLSE. In fact, both the DST and the DNLSE 
contain the same number of constants of motion. 

The analytical solution of eq.(l) in general are not known. However, for nonlinear quan- 
tum dimers which are two sites systems with the nonlinearity either on both the site-energies 
or in one of them can be solved analytically for any arbitrary initial condition. From analyt- 
ical solutions a self tapping transition is found in this model. The trapping of hydrogen ions 
around oxygen atoms in metal hydrides and the energy transport from the absorption center 
to the reaction center in photosynthetic unit have been modeled by the effective quantum 



nonlinear dimer ||l^-p6|. The nonlinear dimer analysis has also been applied to several ex- 
perimental situations, like the neutron scattering of hydrogen atoms trapped at impurity 
sites [|27|J2r 



The self-trapping of a quasi particle at the nonlinear site in a perfect lattice containing 
a single nonlinear impurity has been studied. In the analysis /(| Cm |) has been taken to 
be I Cm The critical value of the nonlinear parameter, x for the self trapping is found 
to increase with increasing the dimension ||29|. It has also been shown that near the self- 



trapping transition the dynamics of the quasiparticle is mostly confined to the few near 
neighbors ^0|- interesting experimental example in this context is the observation that 
trapped hydrogen atoms in metals like Nb move among the sites in the neighborhood of 



impurity atoms such as oxygen. |]15i|27| , ^ Self trapping transition has also been obtained 
in systems where a nonlinear cluster is embedded in a perfect lattice. In fact, two types 
of transitions are obtained, with cluster-trapping transition preceding the self-trapping 
transition . 



Another important feature of this type of nonlinear equations is that these can yield 
stationary localized ( SL ) states and soliton-like solutions. However, the presence of disorder 
or aperiodicity either in site energies or in nonlinearity parameters is needed for this purpose. 



For example, the Ablowitz-Ladik like equation can yield both localized and soliton 

like solution in the presence of disorder in static site energies. Similarly, the DNLSE can 
produce soliton-like solution if any aperiodicity is present in the static site energies PU 



The effect of disorder in the static site energies, in the nonlinear parameters and in both of 
these on the solution of the DNLSE has also been studied by numerical integrations. One 
also finds soliton-like solutions in this case. It has also been shown that the presence of a 
nonlinear impurities in the absence of any disorder in the static site energies can produce 
SL states in one dimension. The first study was made using Greens function approach and 
authors considered /(| Cm |) = | Cm P [0- Furthermore, some cursory discussion on the 



two and more nonlinear impurities has been presented. Later, the one nonlinear impurity 
case has been generalized by taking /(| Cm |) =| Cm T formation of SL states is studied 
in one, two and three dimensions ||^l|-|i3|. 

It is obvious from the above discussion that for a complete understanding of the single 
nonlinear impurity problem the effect of increasing of connections on the formation of SL 



states in the DNLSE needs further consideration. So, we plan to study the formation of SL 
states in Caley tree in the presence of a single nonlinear impurity. Along with /(| Cm |) 
= I Cm \'^, the effect of more generalized nonlinearity impurity ( like rotational nonlinear 
impurity ) needs also consideration. Furthermore, a thorough understanding of the nonlinear 
dimer problem in relation to SL states is required. This is a part requirement to understand 
the effect of disorder in nonlinearity parameters in the solutions of the DNLSE. So, we 
study here in substantial detail the phase diagram of SL states due to the presence of a 
nonlinear dimer in the perfect one dimensional chain. In this connection we also synthesize 
the original Green function approach and the ansatz approach [^pl . Finally, we present 
a logical derivation of DNLSE with /(| Cm |) = | Cm T by generalizing the potential of the 
nonlinear pendulum. To the best of our knowledge no such attempt has been made in this 
direction. 

The organization of the paper is as follows. In the forthcoming section we consider a 
nonlinear impurity of the form — x | C embedded in a Caley tree. We next consider 
in sec. Ill, a rotational nonlinear impurity in the Caley tree. In sec. IV we deal with a 
dimeric nonlinear impurity in an one dimensional system. The impurity is of course, of the 
form —X I C l'^. A similar problem with the rotational nonlinear impurity is considered 
in sec.V. The penultimate section deals with a model derivation of the DNLSE considered 
here. Finally we conclude our paper by highlighting the major features and the importance 
of our results in understanding the behavior of disordered systems made up of nonlinear 
impurities and perfect sites. 



II. A SINGLE NONLINEAR IMPURITY OF THE FORM -x \ C m A CALEY 

TREE 

We consider a Caley tree with the coordination number, Z. Its connectivity is, therefore 
K = Z-1. We put a nonlinear impurity at a site and call it the zeroth site. Furthermore, 
we consider here /o(| Cq |) = | Cq I'^ where a > 0. Now, from this site Z branches will 
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emerge. Again from each of these branches K branches will come out. So, the number of 
points in the n-th generation (or n-th shell) is ZK^^^^\ From any point in the n-th shell a 
quasiparticle can move either to K points in the (n+l)-th shell or to a point in the (n-l)-th 
shell if n> 1. Of course, we are considering only the nearest neighbor hopping with real 
hopping element. So, if the nearest neighbor hopping element is same for all sites, every 
point in a given shell will show the identical time evolution of the probability amplitude 
of the quasiparticle. Thus, the set of equations that describes the time evolution of the 



probability amplitude in the present model is ^4 

.dC„. 



dt 



. dCo 

''~dr 



KCn+i + C, 



n— 1) 



n>l 



-X I Co r Co + ^Ci, n = 



(2) 



(3) 



Furthermore, we have 



y oo 
n=l 



(4) 



Since we are interested here in SL states we write C„(t) = exp(— 'ii?t)a„. After introducing 
this form of C„(t) in eq.(2) and eq.(3) we make the following transformations: (i) a„ = 
(ii) E = E I \fK and (iii) x = x/V^. With this transformations we finally obtain 



+ 0n- 



1; 



and 



= -X 1 00 r 00 + 01, 

/'O r + -r. E I 0" i'= 1 



n > 1 
n = 



(5) 
(6) 

(7) 



(8) 



After some standard algebra for | |> 2 we obtain from eq.(5) and eq.(6) 

0„ _ E-xi^ 1 00 r ^o ^0 

-r ^ C„ol^j - C„+i ol^J 

where G„_o(-^) = \sgn{E)Y^^ Cqo(| E |) rp' and sgniE) denotes the signature of E. We 
also have ri = 1^1-^^^ . So, < r/ < 1. For G°o(| E |) we have 

1 1] 



Go,oi\ E I) 



_4 1 
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(9) 



Eq.(8) for n = yields 

^0 



1 - ^^¥^ + G%{E) = I 00 r G%iE) (10) 



So, the energy of SL states is obtained from solution of eq.(lO). We further note that the 
left hand side of eq.(lO) is independent of the sign of E. On the other hand, to keep the 
sign unchanged in the right hand side of eq.(lO) a simultaneous change of the sign of E and 
X is required. Eq.(lO) when written in terms of t] assumes the form 

--T] = -xK\^o\'' sgn{E) (11) 
V 

Similarly eq.(8) after some simple algebra yields for n > 1 

<Pn = 00 (12) 

Now introducing eq.(12) in eq.(7) we obtain 

- (13) 



Using eq.(13) we eliminate t] from eq.(ll) to obtain the equation for | 0o |- This in turn 
yields 



X 

For K = 1, equation (14) then gives 
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{K+ir-{{K- 


1) + 2 


00 




{K-l + 2 


00 


2~)2 



(14) 



X' 



00 l"^"-'^ (1- I 00 h (15) 

and this equation has already been obtained [JT|. We again note that | 0o P< 1. The energy 
(E) of SL states is 

. ^ . l^ . ^ (2K - (K - 1) \ 60 P)^ , , 

E = -sgnix) iv + -) = -sgn{x)^ ^ ^ ^ ' ^ ' ^ 16 

V V^./(i^+l)2-((i^-l)+|0o|2)2 



Equation (16) can be derived using eq.(13) and eq.(14). We note that for a = 0, eq. (16) in 
conjunction with eq.(14) yields 



E = -sgnix) [-^i^^ + + (17) 



Eq.(17) yields the correct result for a single linear impurity in an one dimensional perfect 
system. We also see from eq.(16) that for K = 1, E = —sgn{x) | X II 0o T^^- Then, for 
(T = 2, I 1 = 1 X |. Inasmuch as the minimum absolute value of the energy of a SL state 
is 2+, the critical value of | x | i.e., \ Xcr \= 2. Above this value of | x 1; two SL states are 



obtained ^ . 



We numerically solve eq.(14) for | <;/)o I as a function of x and a. li \ (po \< 1, we obtain 
an admissible solution. Number of such solutions are the number of SL states in the system. 
We then calculate the energy of the SL states using eq.(16). The phase diagram for SL states 
and energies of these states as a function of x and a are shown in Fig.l and Fig. 2 for K=l, 
2, 3, 10 respectively. We note that for cr = we have the standard Anderson localization 
problem. To obtain a localized eigenstate in a Caley tree of connectivity K, the condition 
I Xcr |> has to be satisfied. This can be obtained from eq.(14) by setting | 0o P= and 
also from equation (17) by setting | E |= 2. This is also obtained in our numerical analysis 
as shown in the Fig.l. For K = 1, there are three regimes. When a < 2, there is always one 
SL state for | x |> 0. For cr > 2, there is a critical value of x (| Xcr |) for each cr ( represented 
by solid curve ) below which no SL state appears but above it two SL states ensue. On the 
other hand, for > 1, we find that for any value of a > 0, there is a | Xcr \ below which no 
SL state appears in contrary to the K = 1 case. But for | x |>| Xcr I, "we get two SL states. 
So, there is a critical line ( dashed line for K = 2, dotted line for K = 3, dot-dashed line 
for K = 10 in Fig.l. ) in the ( x, ) plane which separates the no SL state region and the 
region containing two SL states. If one approaches the critical line from above for any a, 
two distinct SL states merge on the critical line and below that line they disappear. As we 
increase K, the separation between critical lines for any a decreases as shown in Fig.l. So, 
I Xcr I becomes virtually independent of the connectivity for very large K. For the energy 
diagram ( see Fig. 2 ) we have taken cr = 1. We see that for = 1 the energy of the SL 
state increases with x- This implies that localization becomes stronger for larger | x I- 



the other hand for K — 2,3, 10 i.e., for K > 1, one of the two SL states increases with \ x\- 
This corresponds to strong locahzation for large | x I • The energy of the other SL state goes 
towards the appropriate band edge energy. This corresponds to weak locahzation for large 

I X I- 

We now present analytical results in favour of our numerical results. The relevant equa- 
tion for this purpose is eq.(14). We denote right hand side expression of eq.(14) by gK,a{^), 
where A —\ 0o We note that for a > 0, gK,a{Q) = = gK,a{i)- We also have gK,a{A) > 
for A in [0,1]. So, there must be a maximum of this function with respect to A for A in 
[0,1]. For (7 = 0, this maximum occurs at A = 0. If the value of left hand side of eq.(14) lies 
above the maximum, no SL states will be obtained. However, for (4/ | x |^) tying below the 
maximum we shall obtain at least two SL states. Our numerical result suggests that gK,(T{A.) 
has one maximum for A in [0,1]. So, by imposing the maximality condition on gK,cT{A) with 
respect to A i.e., by setting dgK,a{A)/dA — 0, we obtain 

2(7 + 3(t{K - 1)A^ + [{K + if + (X^ + 1)(t\A - aK{K - 1) = (18) 

Note that for cr = 0, A = 0. Furthermore, for K = 1, we get A = ^J{cr - '2)/cr. For K > 1 
and cr > eq.(18) has only one real root, RR given by 

RR 

where, D 
and Dq 

We emphasize that for all value of o" > and K > 1, RR is real and lies between and 1. 
So, the equation for the critical line in the (x, a) plane separating the no state and the two 
states regions is 

I Xcr- 1= ^ y\ ^ (22) 

We now analyze the case of K = 1 along the same line. The relevant equation is eq.(15) in 
which the right hand side will be denoted by ga{A). Inasmuch as < A < 1 the domain of 
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:i-K) ^ 2,^l\a-2){K + lf + 2,'''D- ^^^^ 



2 Q^D 
[^^{K + l)\K - 1) + V3{K + ifDof/^ (20) 



,J{K + 1)2(2 - (7)3 + 21(7{K - 1)2 (21) 



ga{A) is [0, oo] if cr < 2. So, for cr < 2 the line y = constant = always intersect the 

curve (7cr(v4) at a single point. This will happen even if | % | is infinitesimally small. So, for 
cr < 2 we shall always get one and only one SL state if | x |> 0. Precisely this is obtained 



in [^. If o" > 2, 5'o-(0) = = (7ct(1)- Furthermore, gai^) is positive semidefinite for A in [0, 
1]. So, there must be at least one maximum for A in [0, 1]. We also note that for cr = 2, the 
maximum of Qai^) occurs at A = 0. Therefore, for K = 1 and a = 0, applying same analysis 
we obtain | Xcr \= 2 below which there is no SL state and above it there is one SL state. 
Furthermore, it can be shown from eq.(15) by maximizing gcr{A) that | Xcr \ separating the 
no state and two state regions for cr > 2 is 

I I- (^_2)(^-2)/4 y^"^) 



This expression for | Xcr \ has been mentioned without any derivation. 



III. A SINGLE ROTATIONAL NONLINEAR IMPURITY IN A CALEY TREE 

We consider here the case where /(| Cq |) = , !^°L =[f- The origin of this particular 



form of nonlinearity has been discussed in [^]. For this case from the appropriate version 
of eq.(lO) we obtain for | 0o P = A 

+ {K - 1)A^ + PA^ + Q = (24) 
where, P= ^^^^ ^ (25) 

and Q = {K-ir/[x'{^,+K)] (26) 

For K=l, that is for a one dimensional system, we find from eq.(24) 

^X^ - 4 A 

So if I X |> 2, one and only one SL state will be obtained. Furthermore, the energy (E) of 
this SL state is 



4 + A2 

9 



However, K > 1, the analytical calculation of the energy of the SL state is quite prohibitory. 
So, we examine the case for K = 2. Our results are shown in Fig.3. We find again two 
regions. In the lower region we find no SL state. On the other hand in the upper region 
there are two SL states. The line separating these regions contains one SL state. We find 
that for I A |— > , I Xcr oo. This should be expected. Another interesting feature is 
that I Xcr I reduces with increasing | A |. For small value of | A |, the fall in | Xcr \ is quite 
sharp with increasing | A |. 

IV. A NONLINEAR DIMER IMPURITY IN AN ONE DIMENSIONAL 
PERFECT CHAIN: IMPURITY OF -x | C TYPE 

Two consecutive sites occupied by nonlinear impurities are labeled as zeroth site and one- 
th site respectively. We consider here, Xo foil Cq) |) = X | C'q I*^ and Xi /(I Ci \) = x \ Ci T- 
All other Xm are set to zero. All Vn,n+i in eq-(l) are set to unity. Since we are looking for 
SL states, we set C„(t) = 0nexp {—iEt) in the appropriate version of eq.(l). From the set 
of equations thus obtained we solve for 0„ by the standard procedure [^. Consequently, 
we obtain 

^n = G%{E)eo<Po + Gl,iE)e^<P^ (29) 

In equation (29) for the purpose of generality we define eo = — x | 0o T — I 0i T- 

Furthermore, for n = 0, and n = 1 eq.(29) yields 

(Po = Gl,iE)eo(Po + Gl,{E)e^(P^ (30) 
and 01 = G%{E) eo 0o + G° i(^) ei <Pi (31) 

We note that G^^d E |) has been defined by equation (9). After some simple algebra it 
can be shown with the help of eq.(30) and eq.(31) that 

0, = (s^n(E) r/)"-i 01, n >1 (32) 

0_l„,| = (sgniE) 7]p 00, | n | > (33) 
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where < rj < 1 and it is defined in the section before. We further note that the form of 
(pn given by eq.(32) and eq.(33) is independent of the form of eo and ei. Since we must have 
J^'^oo I 0n P= using eq.(32) and eq.(33) we find that 



^0 



P + l0iP=l-^' (34) 
When eo and ei are hnear impurities, from eq.(30) and eq.(31) we obtain 



eoei(eo + ei) ± I 2-eoei I ^(eo-ei)2 + 4 
= 2 (eo - 1) 

where E± denotes the energy of the SL state. We note that for eo = ei > 2, two SL 
states appear above the band. The situation reverses for eo = ei < —2. Furthermore for 
Co = ^ 0, we get E± = ±2. In other words, we get one of the band edge states. These are 
consistent with estabhshed results but the exact form of the solution is not presented 



there. In the present model we have eo = ei = — x for a = 0. On the other hand, since 
I (pQ 1,1 (pi \< 1, for a finite if 00, we get eo = ei — >• 0. General solutions of 

the nonlinear dimer impurity problem are expected to show these asymptotic behaviors. It 
should be noted in this context that for a single nonlinear impurity this continuity has been 



rigorously established 



We now consider the case of our nonlinear dimer. We have three equations, namely, 
eq.(30) and eq.(31) and eq.(34). So, if (po and (pi are fully complex quantities, for the 
problem to be well posed we must have (pi = (pQ. Otherwise, the problem is well posed. 
From eq.(30) and eq.(31) we obtain 

y""-' = I ' ^^^^ 

where, y = | ^ | (37) 



XG%{E) 


00 






-1) 


1 + xG 




<Po 





and X = „ (38) 



It is to be noted that for | (po | = | (pi |, eq.(36) becomes an identity. Consider now the 
case where xGoo{E) > 0. This limit can be obtained either by (i) x > 0, -E > 2 or by 
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(ii) X < 0, £^ < —2. It is easy to see that liY^ 1, eq.(36) is not satisfied. So, the 
solution in this Hmit is | | = | |. Since these quantities can be taken real, we have 
00 = ±01- Consider next the other limit where xGq o(-E) < 0. To achieve this limit we need 
either (i) x < 0, £■ > 2, or (ii) x > 0, -B < —2. Now it can be easily shown that if both 
Do = 1- I xG\o{E) II 00 r< and L>i = 1- | xGI^q{E) \ \ 0i \"< the solution of eq.(36), 
is again | 0o | = | 4>i I- Note further that if F > 1 is a solution of eq.(36), Y < 1 should 
also be a solution of eq.(36) [see for example eq.(35)]. So, any constraint on Dq will imply 
similar constraint on Di and vice- versa. On the other hand we may have | 0o It^I 0i I if 
Dq> Q and consequently also Di > 0. We consider again eq.(30) which for xGq^q{E) < is 

1- I xG%{E) II 00 r=| xGloiE) II 01 r V sgn{E){^) (39) 

We can also consider the equation for 0i (eq.(31)). Both equations will yield the same 
conclusion. It is easy to see that for E > 2, eq.(39) tells that (0i/0o) must be positive. 
On the other hand for E < —2, (0i/0o) should be negative. This result coupled with the 
formula for 0„ (eq.(32) and eq.(33)) suggests that SL states for Dq > 0, either have no nodes 
for £■ > 2 or have N (the number of sites in the lattice) nodes for £■ < —2. We further 
note that band edge states in the perfect system have this characteristics. In case of the 
linear eigenvalue problem, we know that there will be only one SL state in this limit. This 
SL state will evolve from one of the band edge states. Since here we are all dealing with a 
nonlinear eigenvalue problem it is in principle possible to have multiple permissible values 
of 00 and 0i. In other words, we may get a set of states. We shall refer to this a symmetric 
set. But all these SL states will contain same number of nodes. Of course, a subset of these 
states will have 0o = 0i and one state in this subset will survive for cr = 0. 

In the linear dimer impurity problem we also get another SL state. This state has either 
a node {E > 2) through the dimer or the node through the dimer disappears {E < —2). We 
should have the same situation with the nonlinear dimer, but instead of having one such 
state we may have a set of states again with the same number of nodes. Since in the second 
set we must have 0i/0o = —sgn{E) \ (0i/0o) |, for this set to appear we need Dq, Di < 0. 
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00 = 


01 1 



This, in turn implies that the second set of states should always have 
states where (pQ = (pi in the first set together with states coming from this antisymmetric 
set will constitute a major part of the phase diagram for SL states. This phase diagram will 
be discussed here. We also analyze if (po ^ (pi is at all possible in the first set of SL states. 
For I 00 1 = 1 01 I from either eq.(30) or eq.(31) we obtain for | E |> 2, 

V± = 1 II / I 11 (40) 

Ixll0o|'^±l ^ ^ 

For positive sign we need ^ = —sgn{x) = sgn{E) while for the other sign ^ = sgn{x) = 
—sgn{E). Furthermore, the equation for determining | (po \ can be obtained from eq.(34) 
and eq.(40). The required equation for | 0o | is 

(1 - 2 1 00 n (x 1 00 r ±1)' = 1 (41) 

and the energy of the SL state (E) is E = sgn{E){r] + -), < ?/ < 1. 

The phase diagram is obtained by numerical calculation (see Fig. 4). Note that for a = 0, 
we have one state for | x |< 2 and then for | x |> 2 second state appears. This agrees with 
established result On the other hand, we find that for < o" < 2 we find a regime (I) 



with one SL state (coming only from symmetric set). For all values of a in this range there 
is a I XcT I which increases with increasing o", shown in figure by dotted line. Two SL states 
exist along this line. Note further that one of these states coming from the antisymmetric 
set. Above this line we then have a regime (IV) where three SL states are obtained ( one 
from the symmetric set and other two from the antisymmetric set ). 

We now consider the case where a > 2. We have a region (II) where no SL state exists. 



Note that similar situation is also obtained for one nonlinear impurity case . This region 
is bounded by a critical line along which one SL state exists ( shown by solid line in Fig.4). 
This state comes from symmetric set. Above this critical line we have a region (III) of two 
SL states. Furthermore, we never get more than two states from the symmetric set. This 
region of two states continues until we reach the critical line ( shown by dotted line in the 
Fig.4 ) due to the antisymmetric set. Along this line we always have three SL states and 
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above this line we have a region (V) of four SL states. We also mention that o" = 2 is the 
special line. Along this line all transitions are taking place. For example, | Xcr | at cr = 2 
below which no SL state is obtained, is 1. This has been contrasted with the published 
result of I Xcr |~ 2.0 [^]. Another | Xcr \ is found to be exactly 8 which come from the 
antisymmetric set. We also emphasize that cr = 2 line is a line of continuity. The energy of 
states as a function x in regions I and IV of Fig.4 is shown in Fig.5. we have taken a = 1 for 
example to analyze this. In region I the energy of the state always increases with increasing 
I X |. In region IV, the energy of one of the two states belonging to the antisymmetric set 
decreases with increasing | x I ^-^id it is asymptotically reducing to E = 2 { band edge state 
). The energy diagram in regions II, III and V is shown in Fig.6. In this case we have taken 
cr = 3. The energy diagram also shows three distinct regions having no state, two states 
and four states. Two consecutive regions are separated by critical lines where one and three 
states exist respectively. We should note that the energy of one pair of states ( one from 
the symmetric set and the other from the antisymmetric set ) decreases asymptotically to 
E = 2 { band edge state ) with increasing | x I- The energy of other pair increases as | x I 
increases. 

We now provide the analytical support for our numerical results. Since we are considering 
the case | 0o P=| 4>i P here, from eq.(34) we get | 0o P= (1 ~ v'^)/'^- Then, from eq.(40) we 
obtain 



In eq.(42) the upper case sign is for symmetric set while the lower case sign represents the 
antisymmetric set. We consider first the symmetric case. We denote the right hand side of 
eq.(42) by gsi(T,r]) 



and T] lies in [0, 1]. We further note that for a < 2, the domain of gs{<J,ri) is [0, oo]. So, by 
applying our analysis of one nonlinear impurity in one dimensional system we obtain that 




(42) 




(43) 
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for cr < 2, the system will produce one and only one SL state for | Xcr |> 0. In other words, 
in this limit | x 1= 0- Precisely this is seen in our numerical calculation. On the other hand, 
for cr > 2, ^(^((T, 0) = = Qs^cf, 1). Furthermore, gs{<7,r]) is positive semidefinite for r] in [0, 
1]. This, in turn implies that there will be a regime where no SL state will exist. However, 
for I X I exceeding some critical value, | Xcr \, we shall obtain at least two SL states. It is 
quite transparent from the form oi gs{cr,rj) that the function has only one maximum with 
respect to rj in [0, 1] . Our numerical calculations also suggest so. Therefore, there will be a 
no SL state region and a two localizes states region separated by a line in the (x, cr) plane. 
Furthermore, only one SL state will exist along this line. To obtain the equation for this 
line we as usual maximize Qsic, rj) with respect to r). This then yields 

(77^2 - ?7 - 1 = (44) 

and the permissible solution for rj, rjms is 



Vms = 7, (45) 

Now from eq.(43) we obtain for cr > 2 

2^/2 

I Xcr h V (^6) 

9s[cr, r]ms> 

For cr = 2, QsC^, rj) — 2 is the maximum permissible value of gs{2, rf) for r] in [0, 1]. Albeit 
this is not a true maximum of g's(a", rf), the left hand side of eq.(42) must attend this value 
to obtain a SL state. Hence, for cr = 2, | Xcr |= 2. Analytical continuation of eq.(46) also 
yields this result. 

For the antisymmetric case the relevant function to analyze is 

g„((7, 77) =77(l-77)t(l + 77)t-^ (47) 

We note that for cr > 0, g'a(cr, ^) — ^ = Qai^r, !)■ The function is also positive semidefinite 
for Tj in [0, 1]. Hence, there will be a maximum of g(^a, r)) with respect to 7] for rj in [0, 1]. It 
is obvious from the structure of g'a(C) ^) that there is only one maximum. Form these we 
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conclude that there will be two regions, one having no SL state and the other containing 
two SL states. Again the line separating these regions will contain only one SL state. To 
obtain the equation of this line we maximize Qaicr, rj) with respect to rj. We thus obtain 

ar]^ + 7] -1 = (48) 

and the permissible solution of r], rj^a is 



z a 



So, the equation for the Xcr line is 



2W2) 

Xcr 1= T (50) 

ga{(^, Vma) 



We again note that for cr = 2,riam = | which in turn yields gaicTyVam) = |- Hence, from 
equation. (47) we get | Xcr \ = 8. This is obtained from out numerical calculation. For a = 0, 
the maximum permissible value of ga{0,ri) is ^ which occurs at rj = 1. Consequently from 
eq.(50) we obtain | Xcr |= 2. This also agrees with the established result for static dimer 
impurities We further note that the analytical continuation of eq.(50) also produce this 
result. We now note the followings. (1) If cr oo, rims ~ Vma = This implies that two 
critical lines will approach each other as a increases. (2) In contrast to the symmetric case, 
the antisymmetric case will have nonzero | Xcr \ for all cr(> 0) separating the no SL state 
and two SL states regimes. (3) The value of | Xcr I for the antisymmetric case is greater 
than the corresponding value in the symmetric case if a is finite. (4) For the antisymmetric 
case, there is no regime with one SL state except the critical line. This is again in contrast 
to the symmetric case but quite similar to the Caley tree case. 

We now consider the possibility of 00 7^ 0i in the region where Dq, Di > 0. Since E > 2, 
—X > 0, we are taking positive x strictly to avoid complication. We obtain from equation 
for (pn, n & N 

xi\ 00 r - 1 01 n = 4^ (51) 

0100 
1 1 

and 71 = ; ; 1— = ; ; (52) 

Xl0o|'^+J^ xl0i|'^+$ 
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Furthermore, we have eq.(34). We emphasize that not only that 0o and 0i have same sign, 
they are real quantities. For cr = 2, we find after some algebra that 

ri^ -T] + - = (53) 
and P = J liv^Hi? (54) 

Since should be real, we must have ?7 < |, Then from eq.(53) we find | Xcr |= | for this 
case. This is shown by a box in Fig.4. We emphasize that to the best of our knowledge this 
result has not been obtained before. So, along the a = 2 line we have (i) | x |< !> no SL 
state, (ii) 1 <| x 1^ f; one SL state, (iii) | <| x |< 8; two SL state, (iv) | x | = 8, a point of 
three SL states and (v) | x |> 8, four SL states. The case of a = 1, 3 and 4 have also been 
examined. No SL state for (pQ ^ (pi is found. It appears that cr = 2 is a very special case. 
Further analysis for arbitrary a will be presented elsewhere. 

V. A ROTATIONAL NONLINEAR DIMER IMPURITY IN A PERFECT CHAIN 



We consider now the case where fm{\ Cm \) =| Cm P /(yl + (a)^ I 1^) for m = 
and 1 and zero otherwise. In other words we have a nonlinear dimer impurities made up of 



so called rotational nonlinear impurities In this case also 0o = ±01 constitute a set of 
solutions. We focus our attention on this set. More general case will be considered elsewhere. 
We have here two cases, namely the symmetric case, where 0o = 0i and the antisymmetric 
case otherwise. For the symmetric case from the appropriate version of eq.(40) along with 
eq.(34) we obtain 

4A2 



'l + V?[^W-{^-vf] = Gl{ri) (55) 



Since the left hand side of eq.(55) is always positive, we need rj > 1/(1 + A) but at the 
same time we need r] < 1. So, the condition on 77 is 1/(1 + A) < 77 < 1. It is obvious 
that the constraint is met for A > 0. On the other hand, C\{1) = 4A^. Since, G^iv) is 
monotonically increasing function of 1] for rj in [0, 1], 4A^ is the maximum permissible value 
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of the function. So, for a SL state to appear we need | Xcr |= 1- Since, there is only one 
intersection of a line y = = constant with the curve GX{rj) for rj in [0, 1] and | x |> 1; 
we shall always get one and only one SL state from the symmetric case if | x |>| Xcr 1=1- 
For the antisymmetric case from the appropriate version of eq.(40) and eq.(34) we obtain 
4A2 



= (l-77)^[AV-(l + ^r] = GK(77) (56) 



Since G^{r]) should be positive semidefinite in the permissible interval of 77 not exceeding 
[0, 1], we need r]min > -^zi- Since for 77 < 1, we must have A > 2. We note now that 
G~^{1) = = G~^{r]min). So, there must be at least one maximum of G^{rj) with respect to 
77. Imposing the condition of maximality on G~^{r]) with respect to 77 we get 

. - ^' + ' (57) 

It can be easily shown that 77^01; < 1- Furthermore 77^00; 1^ Vmin- Then by setting rjmax — 
Vmin — wc again find that A^r — 2. We also note that G^{r]) has only one allowed 
maximum for r] in [rjmin) !]• Hence, above this maximum we shall not obtain any SL state. 
On the other hand for lying below this maximum, we shall get two SL states. The 
equation of the line separating these two regions is 

I Xcr h / (58) 
\JG/^{r}max) 

We further note that eq.(56) and eq.(57) in the limit A — > 00 yields | Xcr \— 8- This is 

expected. 

In Fig. 7 we have shown graphically the existence or nonexistence of SL states for the 
antisymmetric case with A = 4 for example. The plot of Ga(^) shown as dotted curve. 
Along with that we have superposed the lines y = for three values of x shown in 
figure. The line with | x \— 11-313 touches the maximum point of the dotted curve and 
represents one SL state. So any line with | x |< 11-313 will not touch the dotted curve and 
hence there are no SL states for | x |< 11-313 ( for example the upper line is for | x |= 9.65 
). On the other hand the lines with | x |> 11-313 will touch the line twice. Therefore there 
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will exists two SL states for | x l> 11-313 ( for example the third line with | x 1= 14.61 ). 
Therefore it is very transparent that for this case with A = 4, | Xcr |= 11-313. 

We consider next the full phase diagram obtained numerically- This is shown in Fig. 8. 
For I X |< 1 we have a no state region for any A. For | x |> 1, and A < 2, we have a 
region where one SL state exists ( coming from symmetric set ). For | % |> 1 and A > 2, 
we have two regions separated by critical line (shown by solid line) in the (x, A) plane. 
Above this critical line we have three SL states ( one from symmetric set and other two from 
antisymmetric set ) but below this line we have only one state ( coming from symmetric 
set ). On the critical line only two SL states exists ( one from symmetric set and the other 
one from antisymmetric set ). Energy of SL states as a function of x A = 4 is shown in 
Fig. 9. It shows the typical behavior. 

VI. A MODEL DERIVATION OF THE DNLSE WITH GENERALIZED 

NONLINEARITIES 

To start with we consider a model potential, V{x) given by 

= M '"Tf"'"^ r (59) 



When d = 1, we obtain the potential for the nonlinear pendulum. We also know that such 
a system exhibits two types of motion, namely, libration and rotation. Two regions are 
separated by a critical line, called separatrix |^^. Similar situation also occurs for general d. 



When d = 1/2, the periodic part of the potential is saw-tooth type with kinks smoothened 
out. For d >> 1, the periodic part is more like repeated hard wall potential. These are 
shown in Fig. 10 We further note that for 7 — * 0, V(x) ~ ^ | x p'^. For d = 1, we get the 
usual Hook spring. On the other hand for general d the spring deviates from the Hookian 
behavior. 

We consider next a Hamiltonian, H given by 

. N N 2 

H = -Y.[^n{Xn) I r HCnCl^^ + C^C^+i)] + E [/^ + ^niXn)] = H, + (60) 

^ n=l n=l 
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By defining positions (Qn) and momenta (P„) by Qn = {Cn + C'*)/2 and -P„ = 2(C* — Cn)/2 
respectively, we find that 

Hi = W^n{Xn)iP^ + Ql) + + QnQn+l) (61) 

n n 

So, Hi can be said to define the dynamics of a super oscillator with 2N degrees of freedom. 
Furthermore, the frequency of the oscillator in a given direction is coupled to the motion of 
another oscillator defined by Xn 0. On the other hand, we can replace C„(C*) in Hi by 
the annihilation (creation) operator C„(C^) such that [C.|^,C„,] = ih. C„(C^) then defines 
the annihilation (creation) operator of a particle. Under this transformation it transforms 
to the usual tight binding Hamiltonian of a particle moving on a ring of N sites. Parameters 
Xn then define the displacement of lattice points from their respective equilibrium positions. 
Although we are considering here the local displacement of atoms ( optical phonons ), the 
collective motion of atoms in the lattice ( acoustic phonon ) can be incorporated in this 
formalism by making appropriate modifications of H ||4^. Now from Hamiltons principle 
we find that 

OH 

iCm = 2 = UJjn{Xm) + C'm+l + Cm-1 (62) 

dH pi 

xi = ^ = — (63) 

opi mi 

OH 

and -pi = — = uj'iixi) I Ci p +Vi{xi) (64) 

We now assume that the dynamical evolution of xi occurs in a much longer time scale 
compared to the time scale of Q. So, the dynamics of the other system is determined 
primarily by the adiabatic variation of x„ |^. In the electron lattice interaction language, 
the dynamics of the lattice is assumed to be far far slow compared to the dynamics of the 
electron. Since the adiabatic approximation for xi implies pi = for / G A^, the equation for 
determining xi is then 

cu'iXxi) I Ci |2= -Viixi) (65) 

Consider, for example, Vi{xi) = \kixf and uji{xi) = EiXi. Then from eq.(65) we obtain 
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xi = (Ei/ki) I Ci p. Introduction of this in eq.(62) yields the DNLSE where /;(| Q |) =| Q p 
and xi = Ef/ki. 

In principle can be nonlinear function of x„. Furthermore, the overall dynamics 

of Xfi C3ill be quite complicated and yet can be physically realistic. A few example in this 
regard have already been considered in the literature ETl-H^H^. We consider a situation 



where for n & N Vn{xn) is given by eq.(59) and for uj{x) we assume 



ojix) = — sm{2jx] (66) 

7 V J 

Then from eq.(65) we get 

[1 - 52(^^)2 I C \^]tan^{2^x) - 26 \ C \^ tan{2^x) + S^i^^^) \ C (67) 

where, ^-^-^-^ (68) 

Eq.(67) then yields 



tan{2^x) = ^^[1 + ^1 - y^^^] (69) 

where, ^ = ^ " ' ^ I' (^°) 

Then for d = 1 we have y = 1. This in turn yields tan{2'yx) = 5 | C p. Then after some 
simple algebra we get 



x\c 



|2 



l + (f)^\C 



We consider next the other limit where 7—^0. Then from eq.(69) we obtain 



X 



(2rf-l)|Cp (72) 



k 

This in turn yields uj{x) = — Xo- I C* T where a = 2(2d — 1) and 

X. = -[2^+^ I Eo I (L^d^l^)2.-i] (73) 

Then for d = 1, we get uj{x) = x | C p. Introducing this in eq.(62) we get the desired 
DNLSE. 
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We consider another model where we take u!{x) — 2Eqx and V{x) is given in eq.(59). 
Then again from eq.(65) we get 

Eo\C\^^-kd\x l'^''-^) X (74) 

This in turn yields 

I a; 1 = 1 C 1^ (75) 

kd 

Consequently we have u;{x) — — x(cr') | C \"' where, a' — 2/ {2d — 1) and 

x(a') = 2|Eo|^|^ (76) 
Then inserting the uj{x) in eq.(62) we obtain the desired set of equation. 



VII. CONCLUSION 

The formation of stationary states due to a single nonlinear impurity in a Caley tree and 
a dimeric nonlinear impurity in a perfect ID linear system has been studied here using the 
DNLSE. Two types of nonlinear impurities, namely, /(| C |) =| C where a is arbitrary and 
the rotational impurity are considered. Altogether four cases are studied. Important features 
of these problems are thoroughly discussed in the text. Whenever necessary analytical 
arguments in support of numerical results are provided. Some important aspects of this 
paper are elucidated further below. 

In this paper a very useful synthesis of the usual Greens function approach and the ansatz 
approach is made to gain a better understanding of the problems. Because of synthesis we 
are able to derive for the first time many important results. For example, in the case of 
nonlinear dimer with /(| C |) =| C \'^ , at a — 2 and | x |> | a stationary locahzed state 
in which | 0o 1 7^ I 4>i I is obtained. However, for cr = 1,3 and 4, no such Xcr is found. So, 
it appears that for this nonlinear dimer, cr = 2 is a very special case. Of course, further 
analysis is required and this will be presented elsewhere. 
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For the nonlinear dimer problem, | Xcr \ ~ 2 and 3.6 have been reported in the literature. 
These values of Xcr separate the whole of the a = 2 line into three regions, having no SL 
state, two SL states and three SL states respectively. This has to be contrasted with our 
analysis giving three values of | Xcr \, namely, 1, | and 8. So, we have a new | Xcr \ at 
8 separating two SL states and four SL states regions with three SL states at that point. 
Furthermore, our other two values of | Xcr I are lower than the reported values. 

For the monomer problem | Xcr \ — 8 does not exist. However, one obtains | Xcr \ — 2 
which separates the no SL state and two SL states regions. The analytical reason behind this 
is discussed here. Furthermore, the equation for the critical line separating these regimes 
is analytically derived. In fact, equations describing various critical lines are derived from 
general analytical approach. Our numerical results also agree very well with our analytical 
results. We further note that for the dimer problem this | Xcr | is 1. So, it appears that 
for a cluster of N sites, this | Xcr \ will go as This assertion, however, needs thorough 
investigation. 

Another interesting feature of this problem is the possibility of more localized states than 
the number of impurities in certain situations. This happens due to multiple permissible 
values of the amplitude at the impurity sites. Furthermore, the localization length of one 
set of states increases while this length in the other set of states decreases as the strength 
of the nonlinearity parameter, (x) increases. Hence, in contrast to the static impurity case, 
the system here cam assume more than one configuration. 

Physically, the effect of localized states is manifested in transport properties of the ma- 
terial. However, a system containing a finite number of impurities can only produce finite 
number of such states. So, the probability that these states will exert their influence on 
transport properties of the system is very small, if not unlikely. On the other hand for 
a thorough understanding of properties of disordered system comprising of nonlinear im- 
purities and perfect sites, a critical understanding of these problems is essential. So, the 
importance of the problems presented here should be understood in this context. 
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FIGURES 

FIG. 1. Critical value of nonlinearity (xcr) separating the 'No SL state region' from the 'Two 

SL states region' is plotted as a function of a for Caley trees having K = 1 (solid line), K = 2 

(dashed line), K = 3 (dotted line) and K = 10 (dot-dashed line). There is one nonlinear impurity 

in the Calcy tree of the form -X \ C \"' 

FIG. 2. The energy of SL states shown in Fig. 1 for cr = 1 is plotted as a function of the 

nonlinearity parameter (%) for the Caley tree. The solid curve is for K = 1, the dotted curve is for 

K = 2, the dashed curve is for K = 3 and the dot-dashed curve is for K = 10. 

FIG. 3. The phase diagram for the Caley tree with a single rotational nonlinear impurity 

showing the number of SL states as a function of nonlinearity parameter (x) and the parameter 

A. In this case the connectivity of the Caley tree, K = 2. 

FIG. 4. The phase diagram for a perfect chain with a nonlinear dimer made up of impurities 

of the form -X | C* 1°^- The figure shows regions containing different number of SL states in the 

(x, cr) plane. There are one SL state in region I, no SL state in region II, two SL states in region 

III, three SL states in region IV and four SL states in region V. The solid | Xcr I curve arises from 

the symmetric set and the dotted curve arises from the antisymmetric set. At cr = 2 three values 

of I Xcr I exist. These points are shown by boxes. They are 1, 8/3 and 8 respectively. 

FIG. 5. The energy diagram of SL states as a function of x at a = 1 which corresponds to the 

regions I and IV of Fig.4. The solid curve is due to the symmetric set and the dotted curve is due 

to the antisymmetric set. 

FIG. 6. The energy diagram for SL states as a function of x at cr = 3 which corresponds to 

the regions II, III, V in Fig.4. The solid curve is due to the symmetric set and the dotted curve is 

due to the antisymmetric set. The dashed vertical lines touch the solid and dotted curves at the 
critical values of | x |- 

FIG. 7. The function G^(r/) is plotted as a function of rj (dotted curve). The solid lines 

correspond to the y = for | x I = 9.65, 11.313 and 14.61 respectively. Here A = 4 

FIG. 8. The phase diagram for the perfect chain with a rotational nonlinear dimer impurity 

showing the number of SL states as a function of x and A. 

FIG. 9. The energy for SL states shown in Fig. 8 is plotted as a function of x for A = 4. 
FIG. 10. The potential V{x) = k{^^^^^^Y is plotted as a function of x. The solid curve is 

for A; = 1, 7 = 1, cr = 0.5 and the dotted one is for A; = 1, 7 = 1 and cr = 2. 
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